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Abstract 

The pseudoperturbative shifted - I expansion technique PSLET 
[20-25] is employed to study the spectra of 2D H-donors in magnetic 
field. The Zeeman effect is investigated and results are found to be 
excellent compared to those in ref.[14]. 



1 Introduction 
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The most prominent feature of technology in the last decades was the realiza- 
tion of low dimensional structures. This achievement urged a huge body of 
theoretical and experimental research to study different kinds of interactions 
constrained by the new structures [1-14]. Perhaps, one of the most impor- 
tant and simplest is that of hydrogenic impurity in bulk semiconductors in a 
constant magnetic field. In such systems, the energy scales of both coulombic 
and magnetic potentials are altered because of the dielectric function of the 
medium, e ^ 1, and the effective mass of the electron m* jm ^ 1. In semicon- 
ductors, typical values of the effective mass m* and dielectric constant e make 

771* / 771 

the effective Rydberg R*(— 13.6 — - — eV) about 10 3 — 10 4 times smaller and 
the dimensionless magnetic field 7 = huo c /2R* ( in a fixed magnetic field B 
) about 10 4 — 10 5 times larger than those of the free hydrogenic atom [2,13]. 
As a result, the situation where 7 ~ 1 ( i.e., the intermediate magnetic 
field regime, where neither the Coulombic nore the confining potentials can 
be considered small relative to each other ) is of special physical interest. 
Therefore, many methods had been employed to resolve its significance [14]. 

However, all calculations have shown that the donor states strongly con- 
fined by quantum wells can be correctly described by two dimensional (2D) 
hydrongenlike atoms with properly variational parameters [9-12]. We can thus 
study 2D hydrogenic donors in a magnetic field to understand the behavior 
of donors strongly confined by a quantum wells in different magnetic fields. 
Although the Hamiltonian models in these settings can be separable, one is 
still prevented from obtaining an analytically exact solution. Only for the 
two extreme limits of the magnetic field (7 = and 7 — > 00) the problem can 
be treated exactly [1,3]. Recently, Taut [15,16] has obtained ( conditionally) 
exact solutions for specific values of 7, however no ground-state solutions 
were obtained. 

The adiabatic method [1,3] had been employed in 7 > 1 regime. How- 
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ever, calculations based on this method have usually been restricted to the 
ground state and the first few excited states. Other work has been established 
using variational method. Makado and McGill [4] have based their approach 
on that of Aldich and Greene [2] with different set of basis and they got quite 
good results. Martin et al [11] have used a two-point quasifractional approx- 
imation and found excellent interpolation between the weak and strong mag- 
netic field perturbation expansions. MacDonald and Ritchie [13] have used 
a two-point Pade' approximation. The curves they obtained from different 
Pade' approximation are very different so that no regular pattern appears 
and the results become unreliable [4]. Zhu et al [14] have used a power se- 
ries expansion method. Villalbe and Pino [17] have used a finite-difference 
scheme with a linear mesh of up to 2000 points and failed to provide a good 
estimation of the ground state for the 2D hydrogen atom. Mustafa [18] and 
Quigora et al [19] have used the shifted N-expansion technique and their 
results were in good agreement with those of Martin et al [11]. 

Recently, we have introduced a pseudoperturbative shifted-Z ( I is the an- 
gular momentum quantum number ) expansion technique ( PSLET ) to solve 
Schrodinger equation for states with arbitrary number of nodal zeros. It sim- 
ply consists of using l/l as a pseudoperturbation parameter, where I — I — (3, 
and f3 is a suitable shift. The shift (3 is vital for it removes the poles that 
would emerge, at the lowest orbital states with / = 0, in our proposed expan- 
sion below. Our new analytical, often semi-analytical, methodical proposal 
PSLET has been successfully applied to quasi - relativistic harmonic oscil- 
lator [20], spiked harmonic oscillator [21], anharmonic oscillators [22], two 
dimensional hydrogenic atom in an arbitrary magnetic field [23], two elec- 
trons in quantum parabolic dots [24], truncated coulombic potentials [25], 
...etc [26]. 

Encouraged by its satisfactory performance, we extend, in this paper, 
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PSLET recipe ( in section 2) and introduce its 2D version to treat the problem 
of 2D hydrogenic donors in magnetic fields. In section 3, we discuss the results 
of PSLET. We conclude our work in section 4. 



2 Hamiltonian model and PSLET recipe 



2.1 Hamiltonian model 



The problem of a 2D hydrogenic donor in the presence of a magnetic field 
that is perpendicular to the 2D plane, in the effective mass approximation, 
can be described by the Hamiltonian [14] 



-> he -> jz> we 2 

A ■ A +i — A ■ V 2 

jic ep 



(1) 



where A is the vector potential, u is the electron reduced mass, e is the static 

- d d 

dielectric constant, w is equal to or 1 and V2 represents (i— — h j— ). If we 

ox ay 

choose the cylindrical gauge such that 



A= - B x r 
and the wavefunction \l/(p, if) to be in the form 

Mp) 



*(p,<p) = e 



irrup 



VP 



(2) 



(3) 



then the Schrodinger like equation that describes the system reads 



d 2 m — 1 ( T 2 _ 2w, 

,2 [ \ 



dp 



P- 



P 



u(p) = [E{m) — 7717] u{p) (4) 
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where m is a well-defined magnetic quantum number of the electron in the 
cylindrically symmetric magnetic and 2D Coulomb potentials, E(m) is the 
eigenenergy in units of R* = pe 4 /2h 2 e 2 , the effective Rydberg, 7 = hu c /2R*, 
uj c (= eB/pc) is the cyclotron frequency and p is in units of the effective Bohr 
radius (a* = eh 2 /e 2 p). 

2.2 PSLET recipe 

The 2D version of PSLET starts with shifting the magnetic quantum number 
through / = I — /3, where / = \m\, and (3 is a suitable shift introduced 
technically to avoid the trivial case / = 0, and to be determined below. 
Eq.(4) thus becomes 



d 2 + (j + /?+l/2)(Z + /?-l/2) + ? ■ 



dp 



Q 



u(p)=eu(p) (5) 



~l 2 p 2 an d Q is a scaling factor that 
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where e = E{m) — rwy, V(p) 

is set equal to I 2 at the end of the calculations . 

We employ Taylor's theorem in expanding Eq.(5) about p Q , an arbitrary 
point on the p axis. It is convenient then to transform the coordinates in 
Eq.(5) via x = l l ^ 2 {p — p )/p . Expansions about x — ( i.e., p = p Q ) yield 



+ 



dx 2 



+ V(x(p)) 



u[x) 



Pi 



eu{x), 



where 

V(x(p)) = p 2 J 



1 , V{ Po 



3 + 



Po + Q 

v"( Po )pi 



+ T 1 ' 2 



-2 + 



v\pq)pI 
Q 



X 



2Q 



x 2 + 2/3^(-l) n (n + l)a: n / 



nj—n/2 



n=l 



(6) 
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+ E 



n=3 



:-iy>+ ik + 



'd n V( Po )\ p 2 (p x) 



dp- 



n\Q 



J-(n-2)/2 



i oo 

+ (P 2 - -) E(-!)> + l)x n r {n+2)/2 + 2(3, 



(7) 



?1=0 



and the prime on V(p ) denotes the derivative with respect to p . It is 
convenient to expand e as 



* = E 



(8) 



n=-2 



Equation (6) is exactly of the type of Schrodinger equation for one - dimen- 
sional anharmonic oscillator. This leads to the identification of the leading 
correction in the energy series, Eq.(8), namely 



V{p ) 



E ( ~ 2) = — + 

n >* pV q 



(9) 



Here p Q is chosen to minimize E^ f , i. e. 



dE^ ? d 2 E^ ? 

— - and — > 0, 



dp 



dpi 



(10) 



which in turn gives, with / = y/Q, 



l-P = y/fiV'(p )/2. 



(11) 



Eq.(ll) is an explicit equation in p Q . However, closed form solutions are 
usually hard to be found, if not impossible. Thus we use numerical solu- 
tions to solve for p Q ( hence, the technique is often called semi-analytical 
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). The shifting parameter (3 is determined by choosing lE^ 1 ^ = 0. This 
choice is physically motivated. It requires not only the agreements between 
PSLET eigenvalues and the exact known ones for the harmonic oscillator and 
Coulomb potentials but also between the eigenfunctions as well. Hence 



P = -\(n„ + \)Sl, 



(12) 



where 



n = 2 



3 + 



PoV"(Po) 
V'(Po) 



Equation (7) thus becomes 
V(x(p))=p 2 J 



1 + V(p ) 



pl Q 



+ X>W(:r)r n / 2 , 



n=0 



(13) 



(14) 



where 



(15) 



v (1 \x) 



6Q 



(16) 



and for n > 2 



v (n) (x) 



(-l) n 2 (n + 1) x n + (-1)" (P 2 - l/4)(n - 1) x (n_2) 
+5 n / +2 , (17) 
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(18) 



Equation (6) thus becomes 



d 2 



dx 2 



+ 5>W(x)r n / 2 



n=0 



u(x) = p, 



5>. 

.71=1 



(n-1) 



u(x) 



(19) 



Up to this point, one would conclude that the above procedure is nothing 
but an animation of the eminent shifted large-N expansion ( SLNT ). How- 
ever, because of the limited capabilities of SLNT in handling large -order 
corrections via the standard Rayleigh-Schrodinger perturbation theory, only 
low-order corrections have been reported, sacrificing in effect its preciseness. 
Therefore, one should seek for an alternative and proceed by setting the wave 
functions with any number of nodes as 

u(x) = F npj i(x) exp(U npi i(x)) (20) 
Eq.(19) readily transforms into the following Riccati type equation 

- F:jx) + 2F^ l (x)u' n Jx)] + F nptl (x){ - [U'^jx) + (U' n Jx)) 2 } 



(21) 



+2(3 + -n 2 x 2 + En=i v^(x)l- n / 2 } = p 2 F np j(x) E^l- 



where primes denotes derivatives with respect to xlt is evident that (21) 
admits solutions of the form 



oo n p — 1 

W ( X ) = x n > + E E <l* p i- n/2 

n=0 p=0 



(22) 
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KM = E U^{x)l^ + £ G<${x)T-WI\ (23) 

n=0 n=0 

where 

n+1 

= E^«^ 2m_1 ; A),n,», = 0, (24) 
m=0 

n+1 

= E On,n,n p :r 2m (25) 
m=0 

Substituting equations (22)-(25) into equation (21) implies 

E (i^w (x)z-»/2 + gw (x)/-^ 1 )/ 2 ) -EE (c/^(x)c/(--)(x)r"/ 2 

n=0 P P n=0m=0 P P 

+G£ m >(x)Gk n - ro >(x)H n+2 )/ 2 ^ e « (B) H 2 

p p p p n=0 

n=l n=0 

< ) (# ( * l+1)/2 )]-cw = ° 

(26) 

Obviously, the solution of equation (26) follows from the uniqueness of power 
series representation. Therefore, for a given n p we equate the coefficients of 
the same powers of / and n, respectively. One can then calculate the energy 
eigenvalue and eigenfunctions from the knowledge of C mt n tnp , D mtritnp and 
a^lp in a hierarchical manner. Nevertheless, the procedure just described 
is suitable for a software package such as MAPLE to determine the energy 
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eigenvalue and eigenfunction up to any order of the pseudoperturbation se- 
ries. 

Although the energy series Eq.(8), could appear divergent, or at best, 
asymptotic for small / one can still calculate the eigen series to a very good 
accuracy by performing the sophisticated [N,M] Pade' approximation [25], 



to the energy series, Eq(8). The energy series, Eq(8), is calculated up to 



<V^ 8 by 

e np ,i = P<? + + • • • + <Vr s + 0(1/P), (28) 




(27) 



and with the 



(1//) Pade' approximant it becomes 



[4,4]=P<? + ^4 4 (l/0- 



(29) 



Our method is therefore well prescribed. 



3 Results and Discussion 



For the potential in hand, Eq.(13) reads 



n = 2V2 




(30) 



Eq(ll) in turn yields 




(31) 
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Once Eq.(31) is solved for given values of 7, w, n p and I, the coefficients 
C m ,n,n p , D mjritnp and are obtained in a sequential manner. Then the 
eigenvalues and eigenf unctions are calculated in one batch for each value of 
n p , I, m, 7 and w. 

It worths mentioning that the method is able to reproduce the exact 
know results for both pure Coulombic (7 = 0) and pure magnetic ( w — ) 
interactions. This is evident from Eq.(9). For the pure Coulomb interaction 
( 7 = ), Eq.(9) reads 

P4;3> = -(n, + M + l/2)- 2 , (32) 

while for the pure magnetic interaction ( w — ), Eq.(9) reads 

pEiJ = (2n p + m+\m\ + l)>y (33) 

with identically vanishing higher order corrections. 

Table 1 shows the results of PSLET for the energy spectra of Is, 2p, 2s, 
3d, 3p, 3s and 4d states excluding the Zeeman term iwy, for wide range of 
the dimensionless magnetic field 7. We have plotted the quantum levels of 2D 
donor states without the Zeeman term in Figure 1. It is evident that as 7 
approaches zero, the quantum levels approach those of Eq.(32). Our results 
compare exactly with those obtained by Zhu et al [14]. It is evidently seen in 
Figure 1 that the energy of the quantum levels increases with 7. Increasing 
7 results in a narrower quantum well and hence stronger binding energy. It 
is interesting to note that the energies of the excited states increases more 
rapidly with 7 than those of the ground states and the energy differences 
between these states increase with the applied field. 

To study the effect of the Zeeman term, rrvy, we have obtained the energy 
spectra for small values of the 7. These are displayed in Table 2. Figure 
2 shows the results in table 2 including the Zeeman term. Clearly as the 
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magnetic field is turned on, the states with m/0, i.e., 2p,3d,3p,Ad.., split 
because of their negative and positive m. As 7 increases continuously, the 
energy values of positive and negative m state close to those of lower and 
higher m = states and the levels with negative m cross those of positive m. 
For example, as 7 is turned on the 3d states split into 3d + and 3d- levels. 
As 7 increases, the 3d- state gets closer to the 2s state and crosses the 
2p + state. While the 3d + state gets closer to the 3s state and crosses the 
Ad- state. Clearly 2p_ and 3d- states admit minima due to the competition 
between the Zeeman term and the parabolic quantum well in this weak range 
of 7. 

Figure 3 shows the effect of the Zeeman term at high magnetic field. 
Clearly the energy levels at high magnetic field are accumulated between the 
corresponding m = states and the Landau levels. The order of the levels 
is Is, 2p-,3d- ... for the first Landau level and 2s, 2p +) 3p-, Ad-, .. for the 
second Landau level, and so on. As 7 decreases, the values decreases until 
the hybrid, mentioned above occurs. 

4 Concluding remarks 

Using PSLET, we have obtained the energy spectra of the 2D hydrogenic 
donor in magnetic field. Calculated results have shown the quadratic effect of 
the magnetic field in partially lifting the degeneracy of 2D hydrogenic donor 
states. The linear effect of the magnetic field, i.e., the Zeeman term, has com- 
pletely lifted the degeneracy causing hybrids between different states. The 
order of levels is obtained at high magnetic field. The results obtained here 
show exact agreement with those of Zhu et al [12,14]. The nonperturbative 
nature of PSLET makes it possible to treat the problem irrespective of the 
range of magnetic field. Moreover, we the wavefunctions obtained by PSLET 
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could be manipulated to be used as trail ones in variational approaches. The 
present problem suggests the investigation of excitons in parabolic quantum 
dots as they are known to have same behavior as 2D donors in magnetic 
field. We believe that all aspects of this problem could be revealed through 
PSLET. 
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Figures Captions 

FigllFG 1 ! The 2D donor energy in effective Rydberg units, excluding 
the Zeeman term 7727, is shown as a function of the normalized magnetic field 
7 ranging from to 1 ffor Is, 2p, 2s, 3d, 3p, 3s and Ad states. 

Figure 2: The 2D donor energy in effective Rydberg units, including 
the Zeeman term 777,7, is shown as a function of the normalized magnetic field 
7 for 2s, 2p, 3d, 3p and 4d states. 

Figure 3: Same as figure 1 but including the Zeeman term 7777. 
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